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Reflection and Transmission of Gamma Radiation by 
Barriers: Semianalytic Monte Carlo Calculation j 

Martin J. Berger and John Doggett 

The transport equation for photons which have been Compton-scattered repeatedly with 
a specified sequence of energy losses and deflections was solved analytically, by an "orders of 
scattering" approach, to obtain conditional reflection and transmission probabilities. These 
probabilities were then averaged over all possible intermediate angular and energy configura- 
tions in a Monte Carlo calculation carried out on the NBS automatic computer (SEAC). 
The efficiency of the Monte Carlo calculation was increased by taking advantage of the close 
correlation between the transmission of radiation through thick barriers and the correspond- 
ing penetration problem in an infinite homogeneous medium for which an exact solution was 
available. The reflected and transmitted energy flux (integrated over all spectral energies 
and directions) has been calculated for 0.66-, 1- and 4-Mev radiation incident on water barriers, 
and 1-, 4-, and 10-Mev radiation incident on iron, tin, and lead barriers. Angular and 
energy spectra were obtained for water at 0.66 Mev. The ratio of the scattered energy flux, 
transmitted through a barrier to the scattered flux at equal depth in an infinite medium, was 
found to be constant for distances from the source greater than four mean free paths of the 
source radiation. 



1. Introduction 

The theory of the diffusion and penetration of 
gamma radiation in an infinite homogeneous medium 
is now well developed, complete solutions being 
provided by the moment-method [1] 2 and the closely 
related semiasymptotic Laplace-transform method 
[2]. A satisfactory theory of boundary effects is still 
lacking. The "orders of scattering" approach [3] 
provides formal solutions to boundary problems, 
which can in practice he evaluated only when few 
orders of scattering need be considered. The Monte 
Carlo method [4] is a very effective means for treating 
boundary effects hut is inefficient when applied to the 
radiation that penetrates very far from the source. 
This inefficiency can, in principle, be removed by 
so-called "biased" sampling, but the success of such 
techniques has not yet been demonstrated con- 
clusively. 

This investigation is part of an exploratory Monte 
Carlo program for the survey of boundary effects in 
gamma-ray diffusion. The specific aim is to calcu- 
late the reflection and transmission of radiation by 
plane parallel barriers. The approach is eclectic, 
and a combination of all three of the methods de- 
scribed above is being used. It can be characterized 
as an "orders of scattering" method in which the 
multiple integrals, which ordinarily make the evalua- 
tion of the solution troublesome, are evaluated 
numerically by random sampling. Conversely, one 
could call it a "reinforced" Monte Carlo calculation 
whose efficiency for dealing with deep penetration 
has been increased in two ways: by confining the 
random sampling to the angular and energy variables, 
while treating the space variable analytically, and 
through the use of correlated sampling, the cor- 
related problem being the diffusion of 7-radiation in 
an infinite homogeneous medium. 



1 This work was supported by the Office of Naval Research and the Reactor 
Division of the U. S. Atomic Energy Commission. 

2 Figures in brackets indicate the literature references at the end of this paper. 



The reflection and transmission problems differ 
substantially from the point of view of random 
sampling, because transmission through thick bar- 
riers, being unlikely, tends to require an excessive 
amount of sampling and thus puts a high premium 
on analytical help. Nevertheless, the two problems 
were treated together because they fitted in a single 
computational scheme. 

Two types of questions may be asked in radiation 
diffusion problems. One concerns the estimation of 
the radiation intensity at any one position, cor- 
responding to the observable response of some 
integral detector like an ionization chamber. The 
other type of question is more detailed and concerns 
the spectral composition and directional distribution 
of the radiation at any one point. In the work 
reported here most of the effort has been devoted to 
intensity (energy flux) determinations. Exploratory 
studies of spectral and angular distributions were 
made only for one source energy (0.66 Mev) and one 
barrier material (water) . 

For the deep-penetration problem the attention 
was focused on the comparison of the radiation 
intensity behind a finite barrier of thickness t and 
inside an infinite medium at a depth z (equal to t) 
from the source plane (see fig. l,a,b). No quantita- 
tive indication was previously available on this 
comparison. Notice that practical situations are 
usually intermediate between those schematized in 
figure l,a,b, because there is some matter (possibly 
of different composition) behind the barrier, though 
often at some distance from it (fig. l,c). 

On the problem of relative intensity a rather 
definitive answer has been obtained. In the first 
place, it proved unnecessary to have extended the 
calculation to very great depths because the in- 
tensity ratio no longer varies appreciably beyond 
about four mean free paths of the incident radiation. 
It follows that any future investigation may confine 
itself to moderate depths and may thus utilize less 
elaborate methods. Data on the values of the 
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Figuee 1. Schematization of the barrier reflection-transmission 
problem. 

intensity ratio for various depths are given in table 3 
for four materials (water, iron, tin, and lead) and 
four 7-ray energies (0.66, 1, 4, and 10 Mev). Those 
data may be applied as correction factors to the 
extensive results of calculations for infinite media 
[5]. The ratios differ substantially from unity only 
for low energies and low-Z materials, and also vary 
smoothly, so that adaptation to other circumstances 
ought to be easy. 

In the treatment of the reflection problem, atten- 
tion was directed in the first place to the dependence 
of the reflected intensity on the 7-ray energy, on 
the material of the barrier, and also on the angle of 
incidence. The results for a semi-infinite barrier 
are given in table 2, a. An additional item of 
interest was the dependence of reflection on the 
barrier thickness, that is, on how rapidly the re- 
flection approaches the value for a semi-infinite 
barrier. Some data on this aspect are given in 
table 2,b. 

2. Outline of the Calculation 

2.1. Analytical Frame 

We consider a homogeneous plane-parallel layer 
of material located between the planes 2=0 and 
z=t. Kadiation is incident with specified energy 
and direction on the face z=0. What are the 
characteristics of the radiation reflected through the 
plane 2=0 and transmitted through the plane z=t? 

The state of a photon traversing the barrier is 
described by three variables: the energy E, the 



angle 6 between the direction of motion and the 
2-axis, and the depth z in the barrier. 

Photons can be scattered and absorbed inside 
the barrier. Let ii A (E) equal the probability of 
absorption per unit path length; c(E) the probability 
of scattering per unit path length, and ^(E.6; E' ',0') 
the probability distribution (normalized to unity) 
that a scattering will change the state of the photon 
from (E,B) to (E'fi'). Consider a photon that, as 
the result of n scatterings, has successively assumed 
the energies and directions 
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00, #1, 



> En \ = 

A J 



energy-angle history h n . 



The index n refers to the state of the photon imme- 
diately after its nth scattering; n=0 pertains to the 
state in which the photon enters the barrier. The 
probability distribution of h n is 



Qn= n HEmfirr 
m=0 



E m +\fl m -\-\) 



(1) 



Let P n (Zjt,h n ) be the probability that a photon with 
history h n will cross the plane z (0<z<t) after the 
nth and prior to an (n+l)st scattering. An an- 
alytical expression for this crossing-probability is 
derived in section 3. If a beam of K photons per 
unit area and unit time is incident on the barrier, 
the flux per unit time at depth z is 

F(E,0 9 z,f)=Kibfn fei (dd, . . . (dE n 

n = J J J 

jde n P n QJ(E-E n )5(d-6 n l (2) 

where 8 is Dirac's delta function. If the weight 
factors j n are unity, the flux represents the number 
of photons crossing a unit area of the 2-plane per 
unit time. Most radiation detectors do not measure 
the number flux through a fixed plane; their response 
depends in various ways on the energy and direction 
of the radiation and is represented by the f n 's. A 
few examples are: 



Type 


fn 


Flux 


a 
b 


1 

E n 


Number J Flux (^current) per unit 
Enpre-v area of surface Perpen- 
unergj ^ dicular to the z-axis. 








fFlux per unit area of sur- 








face perpendicular to the 








direction of motion of the 


c 


|sec 6 n \ 


Number 
Energy 


radiation. This is the 


d 


|sec 6 n \E n 


flux seen by an isotropic 








detector. It is propor- 








tional to the photon 








{ density. 



The numerical results of this paper pertain to flux 
type (d). By this choice they become comparable to 
the results of an extensive set of calculations for 
infinite media b}^ the moment-method [5]. 
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2.2. Monte Carlo Estimates 

In the evaluation of expression (2) for the flux, the 
sum over collisons must be cut off for n^>N, and the 
integration over the intermediate energies and angles is 
carried out by the Monte Carlo method. A suitable 
value of the cutoff N depends on the absorbing prop- 
erties of the medium, the depth of penetration, and 
the spectral-energy range over which one wants to 
know the flux. Calculations for an infinite medium 
show that for all practical purposes the energy flux 
vanishes below 30 kev in water, below 50 kev in 
iron, 75 kev in tin, and 100 kev in lead, for source 
energies between 0.5 and 10 Mev. The photon 
energy histories were therefore terminated at these 
limits. Trial calculations indicated that under these 
circumstances, for a barrier thickness up to 16 mean 
free paths (/z £=16), it was sufficient to let N=25 
for water, 12 for iron, 9 for tin, and 6 for lead. 

Energy-angle histories of length N were sampled 
according to the distribution function Q N . The 
detailed equations for doing this are well known. 
They are presented elsewhere, together with the 
adaptations required or desirable for computations 
on an automatic computer [6]. 

Let h nj be an initial portion of length n<N of the 
jth sample history h Nj , and let P nj (z,t) be the corre- 
sponding crossing-probability. If a sample of J 
histories is available, the flux of radiation in a 
specified spectral interval R E and angular interval 
Be can be estimated as follows : 

C C K J N 

dE\ ddF(E,6,z,t)=^J]j: 

J Re J Re J j=ln = 

fnjPnj(z j t)A 1 (E nj )A 2 (6 nj ) } (3) 

where 

\ (W s \f = * > ^ E ni * s m tne interval R E 

ai ^«^ = ^ if En3 is outside the interval R E 

A {R x f = 1 . if nj is in the interval R e 
a^h/^O, if 6 nj is outside the interval R Q . 

If one wants to know the value of the flux for fixed 
values of E and 6 rather than in intervals, the 
sampling procedure must be modified so as to be 
conditional not only on the initial photon state but 
also on the final state E n =E and 6 n =d. For this 
purpose one replaces P n byP*^^.!,^.!; E, 6), 
where 



h*= 



Eq, E u . . . , E n _ 1 ,E 



VQ, 171, ... , <7 w _i, 



We have been content to obtain the flux in energy 
and angular intervals because the main body of 
calculations was done in any case for the total 
intensity integrated over all spectral energies and 
directions. 

In conformity with the usage that has become 
common in shielding work, the results for the 
intensity are presented in the form of buildup 
factors B(z,t), which are equal to the ratio of the 



total to the unscattered energy flux : 



B(z,f)= 



dE\ d9F(E,0,z,t) 
Jo Jo 

Kf P 



(4) 



The Monte Carlo estimates of the buildup factors 
are the mean values 






(5) 



where 



N 



^(2,0=^01 sec ^ | Po)" 1 Z)^»y I sec ^ ni |P ni (2,0. 

71 = 

The standard deviation of the buildup factor 
estimate is J~ 1/2 a(z,t) where 






1/2 



(6) 



W The buildup factors (and their standard devia- 
tions) were estimated from the same set of energy- 
angle histories for three specific situations: 
(i) B(0,t) = energy reflection buildup factor 3 
(ii) B(t,t)= energy transmission buildup factor 
(hi) B (t, oo )= energy buildup factor for the flux at 
depth t within a semi-infinite 
medium. 

2.3. Correlation 

The difference between B(t,t) and B(t,co) is pre- 
dominantly due to photons with histories such that 
in an infinite medium they would cross the plane 
z=t, have their directions reversed as the result of 
one or a few collisions in the vicinity of t, and would 
then make another crossing. Thus B(t,co) is — 
loosely speaking — a function of two factors : a prob- 
ability pi(f) of penetration to depth t, and a proba- 
bility p 2 of being turned around in the vicinity of the 
£-plane. On the other hand, B(t,t) only contains the 
factor j>i(t). The greatest statistical fluctuations in 
a Monte Carlo calculation are found in the estimate 
of pi(t), because only a small part of the sampled 
histories provides information concerning deep pene- 
tration. If the same set of energy-angle histories is 
used for the estimate of B(t,<x>) and B(t,t), it is to be 
expected that the estimate of the ratio of these two 
quantities will have a much smaller statistical dis- 
persion than the estimates of the individual buildup 
factors because the ratio no longer contains the 
troublesome factor P\{t). Unscattered radiation 
should be excluded from the comparison because it 

The energy reflection buildup factor has the' ollowing physical significance. 
Consider an isotropic detector placed in the path of a broad beam of radiation, 
and denote the energy flow through the detector by Fo. Then place a barrier of 
thickness t immediately behind the detector. This increases_the energy flow 
through the detector to F. Then 75(0,0 = F/Fo. The quantity #(0.0—1 is related 
to but not identical with the energy albedo, which as commonly defined pertains 
to an energy current. 
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s not affected by the boundaries. We therefore 
subtract 1 from both buildup factors and introduce 
the ratio 



■g (M) -i . 

5(*,»)-l 



(7) 



A Monte Carlo estimate of k has an estimated 
standard deviation J~ 1/2 a k , where 



Or 



<K*,Q <r(*,oo) 



2p- 



and 



B(t,t)—lB(t,«>)- 



.} 



(1/7) g [5,(M)-5(*,<)] [B,(t, °>)-B(t, »)] 

P(0= <r(M)<K*,°°) 



(8) 



(9) 



is the correlation coefficient of B(t,t) and B(t,oo). 
The two buildup factors have a strong positive cor- 
relation so that p~l and the value of a k /k is small. 
For moderate and deep penetrations, B(t,co) is 
practically identical with the buildup factor for an 
infinite medium which can be obtained by a reliable 
moment-calculation. These known results, together 
with an estimate of the ratio k(f), provide an ac- 
curate indirect method for determining the trans- 
mission buildup factor which requires much smaller 
sample sizes than a direct estimate. Thus the 
Monte Carlo and moment-methods complement 
each other : one provides knowledge of the boundary 
effect, the other knowledge of deep penetration of 
radiation. In section 5 the advantage of this cor- 
relation technique is demonstrated. 

3. Solution of the Transport Equation for the 
Crossing-Probabilities 

The crossing-probabilities obey the following 
integral equations : 



if cos n+ i>O. 



I [dz'P n {z') 

n\ JO 



P„+i(*)= 



c(E n ) 



cos 9 



] -^Uz'P n {z') 



if cos n+1 <O. 



g-[iM(E.. + i)-|-c(E»+i)](z-«')/cos0. + i 

(10a) 

g-[M^(B»+i)+c(E»n))(«-2')/eos«»+i 

(10b) 



Po(z> 



p -[am(#o)+c(Eo)]z/cos6 



(10c) 



These equations express the fact that the probability' 
of passage through the 2-plane after n +1 collisions 
is equal to the probability of passage through the 
plane z' after n collisions, times the probability of 
a collision between z' and z f -\-dz f , times the probabil- 
ity of passage from z' to z without further scattering 



or absorption, the triple product being averaged 
over all values of z f compatible with the boundary 
conditions. 

We shall find it convenient to introduce in place 
of the probabilities P n the related quantities 



G n (z)=\sec6 n \P n (z) 



(11) 



which already incorporate the trigonometrical factor 
required for the type of flux which we want to deter- 
mine. The following abbreviated notation will be 
adopted: 

M« = v*a{E?i) + c(E n ) 
c n =c(E n ) 

rj n =COS 6 n Y (12) 

u n =0 if cos 0„>O 
u n — t if cos 6 n <^0 j 

Equations (10a, b, and c) can then be rewritten in 
the form 



O t 



ViJ%+i 



1 _^? 
6 r (^)=— e vo z 



(13a) 



(13b) 



To solve these equations, assume that G n can be 
expressed in the form 



G n (z)=j:Ale 



-^(z-u m ) 



m = 



(14) 



Substitution of eq (14) into (13a) yields the equation 



G n+1 (z)= 



r n An 



Vn+1 W=0 Vn+l^m VmVn + 1 



i e Vm 



e Vn+1 



—e rim 



- (Z-U m ) 



(15) 



so that the functional form of 6? re+1 is the same as 
that of G n . The expansion coefficients are 



A n m 



1 



Vn+l Mn+1 Mm 



-A n m iorO<m<n (16a) 



4S=-Sit + V 



-U m ) 



(16b) 



We have thus obtained inductively the general 
solution for G n (z). Although the solution is rather 
formal and not readily surveyed, it is convenient for 
numerical evaluation with an automatic computer 
because of its recursive nature. 
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It is possible that rj m ij. n+ i— r) n+ ip m =0, so that singu- 
larities occur in the recursion relations (16a, b). 
But this difficulty is apparent rather than real, 
because upon substitution of the coefficients A^ 1 
into the expression for G n+1 the singularities occur 
pairwise and cancel each other. This can be seen 
by a reexamination of (15). Note that the existence 
of the singularity implies that r] n _i_ir) m ^>0, because 
both ji m and n n+ i are positive. It follows that u m is 
equal to u n+1 . Therefore, the factor multiplying 
A%, in eq (15) is 



(2-M n+1 ) 



C n 6 



As 



Vn+l Vm 



Vm^n+l~~Vn+l^m~ 



(•< 



M n +1 
Vn+l 



V m J 



(2-W„+i) 



— 1 



(17) 



>c n (z—u m )e Vm 



(18) 



If the singularities are removed by a limiting pro- 
cedure, the functional form of G n+i is thus changed. 
For the numerical computation on an automatic 
computer we have found it simplest not to use the 
limiting procedure, but to stick to the unmodified 
formal solution, while eliminating singular cases by 
introducing a cutoff for the quantity a = (/z„ + i/?7 n+1 ) — 
(vmhm) in the denominator of eq (16a). Every 
a<^10" 3 was replaced by 10~ 3 . The error incuired 
by this approximation increases with \z—u m \. But 
even for deep penetration it is quite small. For 
example, when a(z-u m ) = 20, [exp(— 0.001X20) — 1]/ 
10~ 3 =20.201, so that the error is only 1 percent. 

4. Machine Computation 

The calculation was programed for the NBS 
automatic computer (SEAC). It was carried out in 
two stages. First, energy-angle histories (in groups 
of approximately 100) were generated by random 
sampling and stored on magnetic wire or tape. 
These histories were then used as input data for the 
second stage, in which the quantities G n (proportional 
to the crossing-probabilities P n ) were calculated. 
Any desired type of flux or buildup factor can be 
obtained by a suitable linear combination of the 
G n J %, the selection being made by an appropriate 
output code. 

The buildup factor code provides for the computa- 
tion of £(0,0, B(t,t) K and B(t,co), for mt=2 n d 
(n=0,l, . . . , 5), d being a disposable parameter. 
Either all 18 or any subgroup of these buildup factors 
can be computed in an integrated operation. The 
standard deviations a(0,t), a(t,t), and <r(t, °o), as well 
as the correlation coefficient p (t) are also contained 
in the output. 

Most of the computational effort was devoted to 
the buildup factors, but limited use was also made 
of output codes for energy spectra and angular 
distributions. 

The amount of SEAC computing time depends 
on the number N of collisions through which one 



follows the photon histories. As has been mentioned, 
it was found sufficient to let N=25 for water, 12 for 
iron, 9 for tin, and 6 for lead. The amount of 
machine time (seconds per history) required for the 
generation and storage of photon energy-angle his- 
tories was 

U = 0A N. (19) 

For the joint computation of 18 buildup factors, 
the machine time (seconds per history and buildup 
factor) 

* 2 = 0.18 AM-0.0205 N 2 (20) 

for 6 < iV< 25. For the computation of an individual 
buildup factor, the machine time (seconds per his- 
tory) 

2 3 =1.3 N. (21) 

It can be seen that the joint computation of several 
buildup factors was more economical. 

A major difficulty in setting up the 'machine pro- 
gram was the scaling problem. SEAC operates only 
on numbers less than 4 in absolute value. But the 
coefficients A%, in the expansion for G n can occasion- 
ally become extremely large because of sampling 
accidents. This difficulty was aggravated by the 
fact that singularities were not removed analytically 
from the recursion relations (16a, b) but by a 
numerical cutoff procedure. This resulted in the 
occurrence of pairs of coefficients with approximately 
the same large absolute value, but opposite sign, 
whose weighted difference must be calculated with 
great accuracy. This required the use of floating- 
decimal-point routines, which, although they were 
specially adapted to the problem, proved to be 
costly in computation time. 

The occurrence of large terms with alternating 
signs in the series for G n provided a convenient check 
on the code and the computer operation. Clearly 
G n , being a positive multiple of a probability, must 
be a positive quantity. It was found by experience 
that because of the alternating character of the 
series lor G ny almost any sort of code or machine 
trouble quickly led to the occurrence of negative 
G n 's. The machine was instructed to interrupt the 
computation whenever this happened. #j* As many 
G n '& are computed per second, the absence of such 
interruptions for prolonged periods was a reasonable 
guarantee that the operation of thef computer was 
error free. 

5. Results 

5.1. Statistical Efficiency 

Before presenting the main body of our results, 
we shall discuss the statistical accuracy of the Monte 
Carlo estimates in detail for a sample calculation 
based on 300 photon histories, pertaining to 1-Mev 
radiation incident normally on an iron barrier. In 
table 1 are listed, in addition to the estimated 
buildup factors and buildup factor ratios, other 
quantities of statistical interest including fractional 
standard deviations and correlation coefficients. 
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Table 1. Statistical analysis of buildup factor estimates for 1-Mev radiation incident on an iron barrier 



Hot 


5(0, t) 


*(0, t) 
£(0,0-1 


W,t) 


a(t,t) 
B(t,t)-1 


'B(Lco) 


<r(t, °°) 


1-P 


k 


k 


o 


B{U oo ) - 1 


0. 5 

1. 
2.0 

4. 

8.0 

16. 


1. 048 
1. 057 
1. 061 

1. 061 
1. 061 
1. 061 


0. 0586 
.0696 
. 0747 

.0747 
.0747 
. 0747 


1. 40 

1. 71 

2. 43 

4. 07 

7. 80 
17. 8 


0. 61 
. 73 

1.00 

1. 58 

2. 90 
4. 87 


1. 53 

1. 89 

2. 68 

4.45 
8. 60 
20. 


0. 43 

.58 
.89 

1.48 
2.68 
4. 58 


9. 2X10-- 9 

2. 2X10- 2 
6. 5X10-3 

4. 9X10-3 

3. 4X10-3 
1. 85X10-3 


0. 755 

. 798 
. 851 

. 890 
. 895 
. 884 


0. 29 
. 20 
. 15 

. 18 
. 31 
. 41 


0.4 

. 3 

15. 2 

43. 6 
67. 2 

115. 



Now let us see what the large amount of analysis 
introduced into the Monte Carlo calculation has 
accomplished. To estimate the transmission buildup 
factor for a barrier with thickness /z £=16 with a 
standard deviation amounting to 5 percent of the 
buildup factor would require on the order of 10 8 
photon histories in an ordinary random sampling 
calculation. But according to table 1 the semi- 
analytical method requires only 8.5 X10 3 histories. 
Taking into account the estimated increase in 
computation time per history required in the semi- 
analytical calculation, one finds that the effective 
reduction factor is ~ 2,500. This is respectable but 
not quite good enough for a computer with the 
capabilities of' SEAC. But one can do better still 
through the use of correlated sampling, i. e., by 
estimating the buildup factor ratio k(t). It can be 
seen from table 1 that, in accordance with expecta- 
tions, the correlation coefficient p is very close to 
unity, and the fractional standard deviation o~ k /fc 
correspondingly small. A 5 percent standard devia- 
tion now only requires a sample of 67 histories. 
Taking into account the increased computation 
time per history, the reduction of effort is 80-fold. 
The over-all improvement factor is therefore 
2,500X80=200,000. 

Th.eJ.ast column of table 1 lists the factor q = 
[<r(t,t)/B(t,t)YJ[o- k /k] 2 , by which the correlated sam- 
pling technique reduced the required sample size. 
Statistical analyses at other energies and for other 
scattering media yield over-all improvements in 
efficiency of about the same magnitude. 



5.2. Buildup Factors 

Reflection buildup factors for semi-infinite media 
are listed in table 2, a, for diverse conditions. The 
sample sizes were such that the quantities 17(0,t) — 1, 
i. e., the increase of the flux due to reflection, had a 
statistical accuracy of 5 percent. It can be seen 
that energy reflection is appreciable only for low 
energies and low-Z materials. Table 2,b, shows the 
dependence of reflection on the barrier thickness for 
some typical cases. The outstanding feature of these 
results is the rapidity with which the reflection of 
energy assumes its maximum value. A barrier with 
a thickness of one, or at most two mean free paths of 



Table 2. Energy reflection 



a 


. Energy Reflection Buildup Factors for Semi-infinite Barriers 


Energy 


Water 


Iron 


1 
Tin Lead 


00 


0° 


60° 


0° 


60° 


0° 


60° 


0° 


60° 
" 1.042 


Mev. 
0.4 
.66 
1.0 
4.0 

10.0 


1.240 
1.148 
1.081 
1.012 
1.002 


1.278 
1.189 
1.153 
1.035 
1.007 


17661 

1.011 


"L~i42 


" 1.022 
1.005 


"l."686 


T6o9 

1.001 





b. Dependence of Reflection on Barrier Thickness a 


Wate 
Do 

Iron_ 
Do 

Tin. 

Lead 


Material 


Energy 




not 


00 


- 


0.5 


1 


2 


r ... 


Mev. 
0.66 
.66 
1.0 
1.0 
1.0 
1.0 


0° 


0. 65 
.61 
.79 
.89 
.95 
.97 


0.88 
.96 
.93 
.98 
.99 

1.00 


0.99 
1.00 
1.00 
1.00 
1.00 
1.00 




60° 

0° 
60° 
0° 
0° 























a Columns 4, 5, and 6 give the values of fc(0) = [£W)-l]/[#(0,<x>)-i]. 
fc(0) is a measure of the relative amount of reflected radiation from a finite barrier 
and from a semi-infinite medium. 

the source radiation is already equivalent to a semi- 
infinite medium. 

Table 3 contains an extensive list of buildup 
factor ratios k(t), for water at 0.66, 1, and 4 Mev, for 
iron, tin, and lead at 1, 4, and 10 Mev. Again, it 
can be seen that the boundary effect, indicated by 
deviations of kit) from unity, is appreciable only at 
low energies and for low-Z materials. It is signifi- 
cant that the ratios stay constant, within the limits 
of statistical error, as the barrier thickness is in- 
creased from aio£=4 to /zo£ = 16. The physical reason 
for this constancy is the establishment of radiative 
near-equilibria in regions far from the source, which 
is well known from the diffusion theory for infinite 
media. Sufficiently deep in the material, both the 
energy spectrum and angular distribution of the 
radiation flux become very slowly varying functions 
of the penetration depth. Under these circumstances 
it is not surprising that the boundary effect, which 
depends on the characteristics of these distributions, 
likewise remains practically unchanged. 
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Table 3. Buildup factor ratio k(t) = [Bm)-l]/[B(t,^)-l] 
for comparing transmission of radiation through a finite bar- 
rier with penetration to an equal depth in a semi-infinite 
medium 



Materhl 


Energy 


not 


■1.0 


1.0 


2.0 


4.0 


8.0 


16.0 


Water 

Iron 

Tin 


Mev 
f 0.66 

u 

1 10 

{ 10 

I io 


0.601 
.661 

.849 

.790 
.890 
.941 

.889 
.941 
.951 

.939 
.941 

.986 


0.663 
.720 

.885 

.798 
.910 
.959 

.911 
.926 
.960 

.951 
.977 
. 990 


0.713 

.754 
.912 

.851 
.923 
.972 

.924 
.955 
.962 

.969 
.982 
.995 


0.783 
.821 
.920 

.890 
.936 
.974 

.935 
.967 
.973 

.975 
.990 
.992 


0. 785 
.828 
.926 

.895 
.932 
.978 

.938 
.974 
.971 

.979 
.992 
.994 


0.784 
.830 
.933 

.894 
.949 
.977 

.946 
.978 
.969 

.982 
.994 
.995 


Lead 

Estimated 
accuracy 

(%) 


— 


±5.0 


±2.0 


±1.5 


±1. 5 


±2.0 


±2.5 



1 Comparison of a barrier with an infinite medium. 

Table 4. Energy-transmission buildup factors 



Material 


Energy 
Mev 


H 


0.5 


1.0 


2.0 


4.0 


8.0 


16.0 


Water 


( 0.66 

is 

I 10 
I 10 

1 io 


1.49 
1.40 
1.22 

1.40 
1.20 
1.07 

1.29 
1.16 

1.06 

1. 20 
1.11 
1.03 


1.96 
1.80 
1.42 

1.72 
1.36 
1.16 

1.56 
1.31 
1.12 

1.35 
1.23 

1.08 


3. 10 
2.72 
1.83 

2.43 
1.72 
1.35 

2.10 
1.63 
1.26 

1.63 
1.44 
1.17 


5. 99 
5.01 

2.60 

4.07 
2.50 
1.75 

3.15 
2.35 
1.59 

2.09 
1.98 
1.40 


13.3 
10.5 
4.21 

7.80 
4.17 
2.80 

5.31 
4.12 

2.75 

2.87 
3.28 

2.17 


39.4 
25.7 

7.20 

17.8 
7.45 
5.85 

10.2 
9.41 
8.22 

4.24 
7.46 
6.47 


Iron. 


Tin 

Lead 





The results of table 3 pertain to a comparison of 
finite barriers with a semi-infinite medium. For dis- 
tances fiot>4:, a comparison with an infinite medium 
would yield practically identical results. Column 3 
shows the buildup factor ratios for not=l obtained 
by a Monte Carlo calculation for transmission and a 
moment-calculation for an infinite medium. It can 
be seen that these ratios are up to 10 percent lower 
than the corresponding ratios for a semi-infinite 
medium. For /jL t=2 there is a similar but extremely 
small reduction, but our data are insufficient for us 
to quote a value of the reduction. 

Table 4 lists the transmission buildup factors for 
the same range of conditions as in table 3. For 
jj,ot=0.5, 1, and 2, the results were obtained entirely 
by the Monte Carlo calculation. For greater barrier 
thicknesses the values quoted are the result of apply- 
ing the correction factors of table 3 to an extensive 
set of calculations for an infinite medium [5]. The 
infinite-medium results are stated to be accurate to 
within 5 percent. The accuracy ot the derived 
transmission buildup factors is therefore of the order 
of 6 to 7 percent. The direct Monte Carlo results 
for shallow penetration have about the same accuracy. 



5.3. Application of the Monte Carlo Method to an 
Infinite-Medium Problem 

For heavy elements the boundary effect is so small 
that the Monte Carlo method does not add much 
new knowledge. Yet it can be utilized as a con- 
venient and laborsaving tool for the extension of 
available results. As an example of such an appli- 
cation we have considered the penetration of radia- 
tion from a plane monodirectional oblique source in 
an infinite medium. This problem can be handled 
by the moment-method, but requires much more 
work than the corresponding problem for plane 
monodirectional sources with zero obliquity, or 
plane isotropic sources. We considered a 4-Mev 
plane oblique source in lead. A set of 100 energy- 
angle histories were selected that yielded buildup 
factors in close agreement with those of a moment- 
calculation for a monodirectional zero obliquity 
source for penetrations up to p, t=16. Then we per- 
formed a rotation of the angular part of these his- 
tories. Provision for doing this is built into the 
SEAC code for the computation of the G n 'a. After 
rotation, these histories can then be used as input 
data for oblique source problems, and the results can 
be expected to be reliable, since the energy part of 
the histories, known to be properly representative, 




Figure 2. Transmission of radiation from a 4- Mev plane 
monodirectional oblique source through a lead barrier. 
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remained unchanged by the rotation. In this man- 
ner we have made quick calculations for obliquity 
angles 00=45°, 60°, 75°, and 90°. As a final check, 
the fluxes from sources with different obliquities 
were summed so as to obtain the flux from a plane 
isotropic source. Again the agreement with a direct 
moment-calculation was close. Thus a fix on the 
accuracy of the Monte Carlo results was obtained at 
two points, and in spite of the small sample size the 
penetration of radiation to very great depths could 
be calculated with an accuracy which we estimate 
to be 5 to 7 percent. The results of the calcula- 
tion are shown in figure 2 in the form of plots of 
the percent-energy-transmission (^buildup factor 
X exp(— not)) versus not, for different obliquity 
angles. 

5.4. Energy Spectra and Angular Distributions 

In figure 3 we show for normally incident 0.66- 
Mev radiation the energy spectra of the scattered 
flux transmitted through water barriers with thick- 
nesses ju £=4 and not =8. The corresponding energy 
spectra at the same depths in an infinite medium 
are also shown. (The histograms pertain to the 
Monte Carlo results, the curves to a moment- 
calculation.) The shaded areas between corre- 
sponding histograms for transmission and an infinite 
medium are a measure of the boundary effect. It 
can be seen that the presence of boundaries removes 
a large part of the low-energy end of the spectrum, 
but leaves spectral components near the source 
energy unchanged. It is noteworthy that the 
spectra are quite similar at both depths. This lends 
support to our argument that the constancy of the 
buildup factor ratios is connected with the occurrence 
of radiative equilibria. 

Figure 4 pertains to the same conditions as figure 
3, and shows the angular distribution of the scattered 
transmitted energy flux, and of the corresponding 
flux in an infinite medium. The shaded areas 
again indicate the boundary effect. The shaded areas 



E -0.66 Mev 
^t- 4 (WATER) 




E o '0.66 Mev 
Mo f8(WATER) 




Figure 3. Energy spectra of scattered radiation transmitted 
through water barriers of thickness /j.ot = 4 and Hot = 8 } for nor- 
mally incident beams of 0.66- Mev photons. 

The unshaded parts of the histograms pertain to penetration through a finite 
barrier, while the shaded areas indicate the additional radiation that would pene- 
trate to a depth t in an infinitely extended medium. The curves accompanying 
the histograms represent the results of calculations by the moment-method for 
an infinite medium. 



for negative values of cos0 represent photons which 
crossed the boundary 2, 4, 6, • • • times, while the 
areas for positive values of cos0 represent photons 
which crossed the boundary 3, 5, 7, • • • times in an 
infinite medium. The relatively small number 
of photons in the latter group indicates that repeated 
crossings of any plane deep in an infinite medium 
are quite unlikely, so that it should be possible to 
calculate the boundary effect in reasonable approx- 
imation by infinite-medium theory, by setting the 
transmitted flux equal to the infinite-medium 
flux integrated over positive values of cos0. The 
angular distributions are similar to the energy 
spectra in that they do not change their shape as 



~i 1 — 

E =0.66 M ev 
/* t =4 (WATER) 



~l 1 

E =0.66 Mev 

/lot = 8 (WATER) 



1.0 -1.0 

COS 

Figure 4. Angular distributions of scattered radiation trans- 
mitted through water barriers of thickness /jLot — 4- and j*oi = 8, 
for normally incident beams of 0.66-Mev photons. 

The unshaded parts of the histograms pertain to penetration through a finite 
barrier, while the shaded areas indicate the additional radiation that would 
penetrate to a depth t in an infinitely extended medium. The angle B measures 
the direction of the transmitted radiation with respect to the normal to the 
barrier face. Shaded areas for positive values of cos represent photons that 
would cross the boundry plane at least three times in an infinite medium. 
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Figure 5. Energy spectrum and angular distribution of radia- 
tion reflected from a semi-infinite water barrier, for a normally 
incident beam of 0.66-Mev photons. 

The angle 9 indicates the direction between the direction of the reflected photons 
and the normal to the barrier face. 
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the barrier thickness is increased from n t=4: to 
fi t = 8. 

Figure 5 shows the angular and energy distribu- 
tions of reflected radiation due to 0.66-Mev photons 
incident normally on a semi-infinite water barrier. 
Both distributions have two distinct peaks. Their 
locations are such that they appear to imply the 
existence of two preferred modes of reflection, one 
associated with a first scattering through approxi- 
mately 180°, the other with a first scattering through 
approximately 90°. 

6. Discussion 
6.1. Comparison with Experiments 

A recent survey of experiments on the penetration 
and diffusion of y-radiation [5] indicates that few 
experiments have been done for plane source radia- 
tion that are comparable with our calculations. 
Moreover, the uncertainty of the experimental 
results is usually at least of the order of magnitude 
of the boundary effect. It would seem that the 
boundary effect could best be checked by a dif- 
ferential experiment in which one measures the flux 
through a detector, both with and without backing 
material behind the detector. 

The penetration of radiation with energies be- 
tween 0.66 and 2.76 Mev through iron slabs has been 
measured by Beach et al. [7]. Their experimental 
values agree to within a few percent with the buildup 
factors listed in table 4. But the experimental 
results include large corrections for air-scattering, 
so that the evidence is not as conclusive as one might 
wish. 

Kirn, Kennedy, and Wyckoff [8] measured the 
penetration of Co 60 radiation (1.17 and 1.33 Mev) 
through concrete barriers. An attempt has been 
made to match their results by a calculation for an 
aluminum barrier with a concrete-equivalent electron 
density. (The photoelectric absorption cross section 
per electron is very nearly the same for concrete 
as for aluminum, except at extremely low ener- 
gies.) The penetration in an infinite aluminum 
medium was found by interpolation from cal- 
culations of Goldstein and Wilkins [5]. These 
results were then corrected for the boundary 
effect with the use of table 3 of this paper. Finally, 
the predicted buildup factors were increased by 
small amounts (not more than 2%) to account for 
the energy response of the experimental detector. 
One more major adjustment must be made prior 
to a comparison of theory and experiment. The 
experimental cobalt source was 1 cm thick, so that 
there was an appreciable chance for Compton scatter- 
ing and resultant energy degradation of the radiation 
in the source. The degraded radiation is absorbed 
much more strongly than the primary radiation while 
passing through the concrete barrier. When con- 
sidering the attenuation at moderate and great 
depths one should renormalize the source strength so 
as to include only undegraded photons emerging 
from the source. An approximate single-scattering 
calculation leads to the estimate that the effective 



source strength is .82±.()5 times the nominal source 
strength. Table 5 shows the experimental and 
theoretical energy transmission (for normal incidence 
of the Co 60 radiation on the concrete barrier). The 
agreement with the adjusted experimental values is 
good. Although the source strength correction 
constitutes a weakness of the comparison (at least 
insofar as verification of the boundary effect is con- 
cerned), it is significant that by adjusting a single 
normalization constant one can get agreement be- 
tween theory and experiment over a wide range of 
barrier thicknesses. 

Table 5. Energy transmission by a concrete barrier 

Comparison of the experimental transmission found by Kirn et al. [8], and the 
predicted transmission. (Transmission=buildup factor Xexp(— moO-) 





Energy transmission, in percent, 
with concrete barrier thickness of— 


5.0 in. 


10.0 in. 


15.0 in. 


17.5 in. 


Experiment, uncorrected . 


39.0 

46.8 
44.3 


11.0 

13.2 
13.7 


2.80 

3.36 
3.41 


1.40 

1.68 
1.69 


Experiment! corrected for the deg- 
radation of radiation in the source. 
Theory .. .. ._ . .. 





The experimental literature on the reflection of 
gamma radiation from barriers is even scarcer than 
that on transmission. We have learned of some 
unpublished measurements by Kirn, Kennedy, and 
Wyckoff 4 of the reflection of Co 60 radiation from 
effectively semi-infinite concrete and lead barriers. 
(The lead barrier actually was a concrete block 
fronted by }{ in. of lead, but this combination was 
in effect equivalent to a solid lead barrier.) By 
interpolation based on the results in table 2 we 
obtained reflection buildup factors which, according 
to the comparison in table 6, are in reasonably good 
agreement with the experimental values. 

Table 6. Energy-reflection buildup factors for Co 60 radiation 
incident on semi-infinite barriers 

Comparison of the theory with the experiment of Kirn et al. 



Material 


Angle of 
incidence, 0o 


Experiment 


Theory 


Concrete 


deg 

i ° 
I 60 

( ° 
\ 60 


1.063 
1.086 
1.012 


1.065 
1.110 
1. 0095 


Lead _. 




1.041 


1.040 



6.2 Comparison With Other Calculations 

In table 7 a comparison is made for buildup 
factors for shallow and moderate penetration depths, 
pertaining to 0.66 Mev radiation incident on water 
barriers. The comparison includes the results of 
this paper, infinite-medium calculations by the 
moment-method [10] and the results of two other 
Monte Carlo calculations, one of them done by a 
direct stochastic analog method [9], the other 



4 H. O. Wyckoff, private communication. 
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Table 7. Comparison of buildup factors obtained by various 
methods for 0.66-Mev radiation in water medium 



Method 


Transmission, 
B(t,t) 


Penetration in 

semi-infinite 

medium, 

B(t,°°) 


Reflection 
from semi- 
infinite 
medium, 
B(0, oo) 


Monte Carlo (this 
paper) 

Monte Carlo (NBS) 
[6] 

Monte Carlo (NRL) 
[9] .. 


M(/ 


1 


2 


4 


1 


2 


4 


0o=O 


0o= 
60° 


1.96 
1.89 
1.92 


3.10 
3.21 
3.30 


5.99 
5.64 
5.55 


2.48 
2.43 


4.02 
4.07 
3. 75 
3.89 


7.15 
6.98 
6.55 
7.37 


1.148 
1.15 
1.14 


1.189 
1.18 


Spencer-Fano mo- 
ment-method [10] _ . 



based on an analytical integration over a sampled 
collision density [6]. The agreement of the various 
methods is on the whole satisfactory. 

6.3. Comments 

We have demonstrated the constancy of the 
boundary effect for thick barriers, and the resultant 
lack of necessity to do deep-penetration Monte 
Carlo calculations, only for plane geometry. But it 
appears plausible that a similar situation will also 
prevail in other geometries since the constancy of 
the boundary effect appears to be rooted in the 
occurrence of equilibrium spectra far from a radiation 
source. 

The advantage of incorporating knowledge ob- 
tained from other sources into a Monte Carlo calcu- 
lation through the technique of correlated sampling 
has been clearly demonstrated. In fact, it was only 
this technique which made it possible for us to survey 
the boundary effect for a variety of conditions with- 
out excessive computation. 

The analytical treatment of the space variable 
likewise proved successful, but we feel that it may 



have been too much of a good thing. From the stand- 
point of computing efficiency, there is an optimum 
amount of analysis that should be put into a Monte 
Carlo calculation, which depends on the nature of 
the problem and the capabilities of the computing 
machine. It is questionable whether we operated 
near this optimum. If one is interested in the 
radiation intensit}^, it is not unreasonable to take, 
as we have done, relatively small samples and to 
analyze them exhaustively. But if energy spectra 
and angular distributions are desired, it is probably 
preferable to use larger and more representative 
samples of energy-angle histories and to analyze 
them less thoroughly. One possible compromise 
would be to confine the analytic treatment of the 
space variable to a limited number of collisions, and 
then to revert to a stochastic calculation. 



The authors thank Anne Futterman for aid and 
advice in programing the problem for SEAC. 
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